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In conventional micromagnetism magnetic domain configurations are calculated based on a con- 
tinuum theory for the magnetization which is assumed to be of constant length in time and space. 
Dynamics is usually described with the Landau-Lifshitz-Gilbert (LLG) equation the stochastic vari- 
ant of which includes finite temperatures. Using simulation techniques with atomistic resolution 
we show that this conventional micromagnetic approach fails for higher temperatures since we find 
two effects which cannot be described in terms of the LLG equation: i) an enhanced damping when 
approaching the Curie temperature and, ii) a magnetization magnitude that is not constant in time. 
We show, however, that both of these effects are naturally described by the Landau-Lifshitz-Bloch 
equation which links the LLG equation with the theory of critical phenomena and turns out to be 
a more realistic equation for magnetization dynamics at elevated temperatures. 

PACS numbers: 75.10.Hk, 75.40.Mg, 75.75.+a 
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An increasing amount of research is focusing on the 
dynamic behavior of ferromagnetic materials at elevated 
temperatures. The motivations for this are manifold. A 
major imperative is the understanding of pulsed laser 
experiments on thin film samples, for example the all 
optical FMR experiments of Van Kampen et. al. [l|, 
and the higher laser power experiments of Beaurepaire 
et. al., [21 who demonstrated complete demagnetization 
on a timescale of picoseconds. One of the main issues of 
the high-temperature magnetization dynamics is the rate 
of the magnetization relaxation due to different processes 
involving magnon, phonon and electron interactions that 
contribute to thermal spin disordering. 

The basis of most of theoretical investigations of ther- 
mal magnetization dynamics is a micromagnetic ap- 
proach which considers the magnetization of a small par- 
ticle or a discrete magnetic nanoelement as a vector of a 
fixed length (referred to here as a macro-spin) with the 
phcnomenological Landau-Lifshitz-Gilbert (LLG) equa- 
tion of motion augmented by a noise term J3]. However, 
contrary to atomic spins, there is no reason to assume 
a fixed magnetization length for nanoelements at non- 
zero temperature. For instance, the latter can decrease 
in time upon heating by a laser pulse. Hence, from the 
point of view of modeling of magnetization dynamics, 
there is a general need for further development of the 
micromagnetic theory in terms of its ability to deal with 
elevated temperatures. 

Within this context we note the failure of micromag- 
netics in general to deal with the high frequency spin- 
waves which give rise to the variation of magnetization 
with temperature. It has been suggested to treat this 
problem using scaling approaches [J, |5( . A similar prob- 
lem arises in multi-scale modeling (with atomistic and 
micromagnetic discretizations to treat, for example, in- 
terfaces p| ) which can not correctly describe the transfer 



of high energy spin-waves from atomistic into the micro- 
magnetic region. An alternative approach is the coarse 
graining model of Dobrovitksi et. al. [TJ, which has the 
advantage of being able to link the length-scales but has 
been developed for simple systems only. 

Some understanding of the pulsed laser experiments 
could indeed be obtained in terms of a micromagnetic ap- 
proach taking into account, in an empirical way, the tem- 
perature variation of the intrinsic parameters, particu- 
larly the saturation magnetization M^ and the anisotropy 
energy density K. Lyberatos and Guslienko jj] have 
used this macro-spin model to investigate the response of 
nanoparticles during the Heat Assisted Magnetic Record- 
ing (HAMR) process. The validity of the macro-spin 
approach including the thermal variation of model pa- 
rameters has further been investigated in Ref. [9j using 
an atomistic approach. This work demonstrates that, al- 
though the macro-spin model works well for temperatures 
far below the Curie temperature Tc, longitudinal fluctua- 
tions of the magnetization become important at elevated 
temperatures, which cannot be treated within the macro- 
spin model of the corresponding LLG equation of motion. 
The use of a macro-spin of fixed length places the same 
physical constraint on micromagnetics at temperatures 
close to Tc. Clearly, some approach to macro-spin dy- 
namics beyond the LLG equation is needed. 

A semi-phenomenological equation of motion for 
macro-spins allowing for longitudinal relaxation has been 
derived in Ref. [iflj within the mean-field approxima- 
tion (MFA) from the classical Fokker-Plank equation for 
individual spins interacting with the environment. This 
"Landau-Lifshitz-Bloch (LLB) equation" has been shown 
to be able to describe linear domain walls, a domain wall 
type with non-constant magnetization length. The valid- 
ity of these results has been confirmed by measurements 
of the domain wall mobility in YIG crystals close to Tc 



[ll| and by recent atomistic simulations jl2| . 

In this letter we explore high-temperature dynamic 
properties using atomistic modeling. These simulations 
are still based on the LLG equation on the atomic level 
and, hence, do still not provide a microscopic description 
of the damping itself. Nevertheless they do include ther- 
mal degrees of freedom microscopically and encapsulate 
important phenomena associated with relaxation, includ- 
ing the thermodynamics of the phase transition and both, 
longitudinal and transverse macroscopic relaxation. We 
find an enhanced transverse relaxation when approaching 
the Curie temperature from below and a magnetization 
magnitude which is not constant in time. Both of these 
phenomena cannot be understood in terms of conven- 
tional micromagnetism but, comparing these predictions 
with a macro-spin model based on the LLB equation, 
we conclude that here these phenomena are indeed well 
described by the LLB equation. 

For our atomistic simulations we use a model in which 
the dynamic behavior of classical spins |si| = 1 on lattice 
sites i with magnetic moment jiq is treated at the atomic 
level with the Langevin form of the LLG equation 



Si = -7[si X Hi] - ja[si X [s^ x H^ 



(1) 



where 7 is the gyromagnetic ratio, and a is the damp- 
ing parameter, a = 0.1 in our simulations. The total 
field Hi contains nearest-neighbor Heisenberg exchange 
(exchange constant J) and Zeeman contributions and 
it is augmented by a white- noise field (i{t) with the 
correlator {Cif,it)QAt')) = ^^S,j5^^ 5(1 - t'), where 
fi,v = X, y, z. For simplicity, the dipolar interaction is 
neglected as well as any crystalline anisotropy. A cu- 
bic lattice with periodic boundary conditions and system 
sizes of 48^ has been considered. In the calculations we 
first establish thermal equilibrium for a given tempera- 
ture starting with all magnetic moments parallel to the z 
axis and applying a field Hz — 0.05 J/ /LtQ. Then, to eval- 
uate the transverse relaxation, all spins were simultane- 
ously rotated by an angle of 30°. We have calculated the 
average spin polarization m = (1/N) J^i i^i) per lattice 
site which is proportional to the experimentally observed 
magnetization M. 

Fig.^ shows one transverse magnetization component 
as a function of time for different temperatures. The 
magnetization is normalized to its initial value and the 
data show clearly a faster relaxation for higher tempera- 
tures. Note that in our simulation even above the Curie 
temperature Tj, there is still a finite magnetization due 
to finite-size effects and the fact that the simulations are 
conducted within an external field. Fitting the curves 
to an expression mx{t) ~ cos(i/rp) exp(i/Tj_) shows a 
perpendicular relaxation time t± which increases with 
temperature, deviating from its zero temperature limit 
l/{ajHz)- Fig. ^3 presents the change of the absolute 
magnetization value as a function of time for a similar 
simulation but with a large angle of 135° . Note that the 
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FIG. 1: Relaxation of the magnetization for different temper- 
atures using the atomistic modeling: a) normalized perpen- 
dicular component (30° excitation) ; b) absolute value of the 
magnetization m = |m| (135° excitation). 



magnetization magnitude shows a dip during the relax- 
ation process which is well below its equilibrium value. 
A dynamic response of this type cannot be described in 
terms of the macro-spin LLG equation which conserves 
the absolute value of the magnetization, but is consistent 
with the LLB equation as will be discussed below. 

It is interesting to note that our atomistic model with 
a constant microscopic damping parameter exhibits an 
increase in the effective macroscopic damping as ob- 
served experimentally [l3j . We believe that this is due to 
magnon-magnon scattering processes which give rise to 
the initial decrease of the magnetization as the energy is 
transferred from the k — (precessional) mode to higher 
order modes. This results in the enhanced transverse 
damping and in the dip of the magnetization, followed 
by the recovery to its equilibrium value as the spin waves 
decay. 

Furthermore, we investigate the longitudinal relax- 
ation time r|| from the initial relaxation of the fully or- 
dered system to thermal equilibrium. The relaxation of 
the magnetization to equilibrium is found to be approx- 
imately exponential on longer time scales which defines 
the characteristic time r||. Fig. |21 shows the variation of 
the longitudinal relaxation time with temperature. The 
rapid increase close to Tc is known as critical slowing 
down 14], a general effect characterizing second order 
phase transitions. Also shown in Fig. [3 is the perpendic- 
ular relaxation time tj_ determined as described above. 
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FIG. 2: Temperature dependence of longitudinal and trans- 
verse relaxation times from the atomistic modeling and the 
LLB equation, calculated as inverse rates given by Eq. ©. 



Approaching the Curie temperature the perpendicular 
relaxation time t± breaks down. 

As we have demonstrated so far, the atomistic model 
shows important physical aspects of the behavior of 
nanoscale magnetic systems, including a temperature 
dependence of the effective damping, longitudinal fluc- 
tuations and critical slowing down. Next, we demon- 
strate that these effects can be described alternatively 
by macro-spin magnetization dynamics in terms of the 
Landau-Lifshitz-Bloch equation of motion [lO| . This pro- 
vides not only a deeper understanding of the phenomena 
but it also suggests that the LLB equation is more suit- 
able than the LLG equation for finite temperature mi- 
cromagnetics. 

The LLB equation following from Eq. (Q in the spa- 
tially homogeneous case can be written in the form 

(m • Hcff)m 



rfi — — 7[ni X Hoff] + jan 



7"-L 



[m X [m X Hcff]] 



(2) 



where m = (s) is the spin polarization and a\^ and a±_ 
are dimensionless longitudinal and transverse damping 
parameters given by 
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(3) 



for T < T^^^ and the same with a± => ay for T > 
rpMFA^ where T^^-'^ is the mean-field Curie temperature. 
Here, a is the same damping parameter that enters Eq. 
(^ . The effective field Hoff is assumed to be much weaker 
than the exchange interaction and it is given by 
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(4) 



Here H and H^ are applied and anisotropy fields and rrie 
is the zero-field equilibrium spin polarization in the MFA 
that satisfies the Curie- Weiss equation 



m = B [f3{mJo + ^aH)] 



(5) 



with F = and e = 1 - T/T^^^. B is the Langevin 
function, P — 1/ (kBT) , and Jq the zero Fourier com- 
ponent of the exchange interaction related to T^ 



MFA 



fcsT^MFA ^ j^/3_ jj^ Eq EJ XII = dm{H,T)/dH is the 
longitudinal susceptibility at zero field that can be ob- 
tained from Eq. jSJ. The anisotropy field Ha due to 
the uniaxial anisotropy is related to the zero-field trans- 
verse susceptibility x± ^s Ha — {m-xex + "m-y^y) /x±lld|- 
The equilibrium solution of the LLB equation satisfies 
m X Hoff = and m • Heff = 0. For T < T^'^^ the lon- 
gitudinal susceptibility x\\ becomes very small in which 
case it can be shown that m = me- This means that the 
longitudinal relaxation vanishes and Eq. (0) reduces to 
the standard LLG equation with a± = a. 

In the damping parameters a\\ and a± of Eq. © a 
is non-critical at T^ . Its temperature dependence can- 
not be established within our semi-phenomenological ap- 
proach, so we assume it to be a constant, for the sake 
of comparison with the results of our atomistic simula- 
tions. The LLB equation also can be written in terms 
of the vector n = m/nie UM ■ This form provides a link 
to the micromagnetic anisotropy constants but becomes 
inconvenient above Tc where me disappears. 

In order to effect a comparison we analyse the relax- 
ation rates derived from the LLB equation. Firstly we 
note from Eq. |Sla linear increase of a|| with T, whereas 
the behavior of a± is non-monotonic, changing from a 
linear decrease below T^^-^ to a linear increase above 
jiMFA^ However, it is important to note that ay and 
a±_ are non-critical for all finite temperatures, and that 
the variation of a±_ is weak. With this background, we 
now consider the relaxation rates from the linearized LLB 
equation which have the form 
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where x\\{H^T) is the longitudinal susceptibility at 
nonzero field that follows from Eq. (O or simply from 
m • Hoff = 0, in our approximation. 

The longitudinal relaxation rate is, in general, very fast 
as Fy ~ Jq. Since x\\{H,T) is large near T^^^, T\\ shows 
critical slowing down which is a result of the critical be- 
havior of Xji (i/,T) rather than the variation of a\\. The 
transverse susceptibility for the isotropic model is sim- 
ply given by x±(i/, T) = m{H, T)/H so that Fj, ~ F is 
much smaller than Fy below T^^^. However, it increases 
with temperature, as was observed in the atomistic mod- 
eling presented above and its critically behavior close to 
Tc is Fj^ ^ l/m{H,T). For temperatures below Tc a 
corresponding behavior was found for the line widths of 
FMR experiments [IJJ . 
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FIG. 3: Relaxation of the magnetization for different temper- 
atures as in Fig. but using the macro-spin LLB modeling. 



At T = T: 
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the rates are given by 
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1 /s 
where toj;^ = [(5/3) {^oH/Jq)] is the induced magne- 
tization at T^MFA^ Above T^^^ both rates merge: 
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Finally, in the presence of uniaxial anisotropy F j_ is given 
by Eq. ® with l/xj.(iJ, T) = H/m{H, T) + l/x±, where 
X_L is only weakly temperature dependent within mean- 
field theory below T^'^^. 

To compare the LLB results with the predictions of 
the atomistic model, Fig. |21 includes the inverse relax- 
ation rates calculated using Eq. IjHJl with rescaled tem- 
perature to fit the exact value fcsTc = 1.44J for a simple 
cubic lattice. The agreement between Eq. © and the 
numerical results is remarkable given the MFA used in 
the derivation of Eq. © . 

Also, we have integrated numerically Eq. Q for a 
macro-spin to give the time evolution of the magnetisa- 
tion components for comparison with the numerical re- 
sults of Fig. ^ The results are presented in Fig. 13 Com- 
parison with Fig. n shows that the LLB equation repro- 
duces essential physical processes which govern the mag- 
netization dynamics at elevated temperatures and thus it 
can be used as an alternative to micromagnetics in this 
region. However, this comparison could still be improved 
if one evaluates the macro-spin parameters directly from 



an atomistic simulation. Furthermore, if the LLB equa- 
tion is to be used as an alternative to micromagnetics, 
the corresponding parameters could as well be extracted 
from experiment. 

In conclusion, performing atomistic simulations of 
thermal magnetization dynamics we observe an increase 
of the macroscopic transverse damping approaching the 
Curie temperature. This increase is determined by the 
thermal dispersion of magnetization and would exist in- 
dependently from any other possible thermal dependence 
of internal damping mechanisms such as phonon-magnon 
coupling. This effect explains the broadening of the reso- 
nance line width in classical FMR experiments [iJl . Fur- 
thermore, the magnetization vector turns out not to be 
constant in length. Instead during relaxation one can 
observe a dip of the magnetization which is more pro- 
nounced when approaching the Curie temperature. Fi- 
nally, the magnetization dynamics has important con- 
tributions from longitudinal relaxation. This relaxation 
shows critical slowing down at temperatures close to Tc- 
Importantly, the observed dynamics is in agreement with 
the dynamics of a macro-spin described by the Landau- 
Lifshitz-Bloch equation which contains both longitudinal 
and transverse relaxation. This equation could serve in 
future as a basis for an improved micromagnetics at ele- 
vated temperature. 
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